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We present results from our studies of radio emission from selected Ultraluminous X-ray (ULX) sources, using archival 
Giant Metrewave Radio Telescope (GMRT) data and new European VLBI Network (EVN) observations. The GMRT data 
are used to find possible faint radio emission from ULX sources located in late-type galaxies in the Chandra Deep Fields. 
No detections are found at 235, 325 and 610 MHz, and upper limits on the radio flux densities at these frequencies are 
given. The EVN observations target milliarcsecond-scale structures in three ULXs with known radio counterparts (N4449- 
XI, N4088-X1, and N4861-X2). We confirm an earlier identification of the ULX N4449-X1 with a supernova remnant and 
obtain the most accurate estimates of its size and age. We detect compact radio emission for the ULX N4088-X1, which 
could harbour an intermediate mass black hole (IMBH) of 10 J Mq accreting at a sub-Eddington rate. We detect a compact 
radio component in the ULX N4861-X2, with a brightness temperature > 10 6 K and an indication for possible extended 
emission. If the extended structure is confirmed, this ULX could be an HII region with a diameter of 8.6 pc and surface 
brightness temperature > 10 s K. The compact radio emission may be coming from a ~ 10 Mq black hole accreting at 
0.1AW 
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1 Introduction 

Several scenarios have been proposed to explain the high 
luminosities (Lx > 10 39 erg/s) of Ultraluminous X-ray 
sources (ULXs), but none of them is able to reveal the phys- 
ical nature of all ULXs. If ULXs are powered by accretion 
at the Eddington rate, this would imply accreting compact 
objects of masses 10 2 -10 5 Mq. Such intermediate compact 
objects can only be black holes (Colbert & Mushotzky 1999) 
and they would be the missing link between stellar mass 
black holes and supermassive black holes in the nuclei of 
galaxies. These Intermediate Mass Black Holes (LMBHs) 
could form from the death of very massive and hot stars 
or from multiple stellar interactions in dense stellar clusters 
(Portegies Zwart 2003). It has also been suggested that ULX 
objects may harbour secondary nuclear black holes in post- 
merger galaxies (Lobanov 120081 ), with masses in excess of 
10 5 Mq and accreting at sub-Eddington rates. Alternatively, 
ULXs could be neutron stars or stellar mass BHs apparently 
radiating at super-Eddington luminosities (Begelman 2002). 

Radio observations of ULXs bear an excellent poten- 
tial for uncovering the nature of these objects, by detecting 
and possibly resolving their compact radio emission, mea- 
suring its brightness temperature and spectral properties, 
and assessing the physical mechanism for its production. 
Few ULXs have been studied in the radio domain (Kaaret 
et al. 2003; Kording et al. [2005 ) and a small sample of 
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ULXs has been cross-identified in the existing radio cata- 
logs (Sanchez-Sutil et al. 120061) . 

An increase of the number of radio detections and sub- 
sequent Very Long Baseline Interferometry (VLBI) studies 
of detected radio counterparts could potentially help to clar- 
ify the nature of ULX sources. With this aim, we 1) ana- 
lyze archive images of the Chandra Deep Fields taken with 
the Giant Meterwave Radio Telescope (GMRT) looking for 
faint radio counterparts of the ULX sources located in this 
fields; and 2) initiate an European VLBI Network (EVN) 
program to detect and study milliarcsecond-scale emission 
in ULX objects with known radio counterparts. In Section 2, 
we present the two samples of ULX objects studied. The 
observations and data reduction are explained in Section 3. 
The results obtained are shown in Section 4, leading to a 
discussion and final summary presented in Section 5. 

Throughout this paper we assume a A cold dark matter 
(CDM) cosmology with parameters Hq — 73 km s _1 Mpc 
fl A =0.73 and Q m = 0.27. 

2 Observing targets 

We use archival GMRT data to search for radio counter- 
parts of 24 ULX objects identified in the Chandra Deep 
Field North (CDFN), Chandra Deep Field South (CDFS), 
and Extended CDFS (Lehmer et al. 120061 . All these ULXs 
have luminosities Lx > 10 39 erg/s in the 0.5-2.0 keV band, 
and are located in optically bright irregular and late-spiral 
galaxies. Ten of the 24 X-ray sources appear to be coinci- 
dent with optical knots of emission, with optical properties 
that are consistent with those of giant HII regions in the 
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local universe, suggesting that these ULX sources trace dis- 
tant star formation (Lehmer et al. 2006). 

The objects targeted in our EVN observations are se- 
lected from a sample of 1 1 ULX objects with radio coun- 
terparts (Sanchez-Sutil et al. l2"0"06j l in the VLA FIRST cat- 
alog (Becker et al. 19951). We select three ULX which are 
brighter than 1 mJy and clearly away from the nuclear re- 
gion in their respective host galaxies. The first target, N4449- 
X4, is identified as the most luminous and distant member 
of the class of oxygen-rich Supernova Remnants (SNRs) 
(Blair et al. |19831 l and classified as an ULX source by Liu 
& Bregman 120051 A detailed study of this source is de- 
scribed in Mezcua & Lobanov (in prep.). The second target, 
N4088-X1, is an ULX located at a distance of 13.0 Mpc in 
the asymmetric spiral galaxy NGC 4088. This ULX is lo- 
cated within the extended emission of a spiral arm and is 
coincident with a conspicous maximum of radio emission 
of 1.87 mJy (at 1.4 GHz) with an offset of 3.62" to the X- 
ray peak (Sanchez-Sutil et al. 2006). The ULX has a X-ray 
luminosity in the 0.3-8.0 keV band of 5.86 x 10 39 erg/s (Liu 
& Bregman 120051 who don't rule out the possibility of it 
being an HII region). The third target, N4861-X2, is located 
in the spiral galaxy NGC 4861 and has an X-ray luminosity 
of 8.4 x 10 39 erg/s (Liu & Bregman I2005I ). Its radio coun- 
terpart is offset from the X-ray position by 1 .97". This ULX 
has been suggested to coincide with an HII region powered 
by massive early OB type stars (Pakull & Mirioni [20021) . 

3 Observations and data reduction 

We use archival GMRT data of the Hubble Deep Field North 
(overlapping with our CDFN region of interest) at 235 MHz 
(experiments 1 NIK04 & 1 1 TMA0 1 ) observed on 200 1 De- 
cember 17, 2002 January 4 and 2002 January 26, and of 
the CDFS at 325 & 610 MHz (experiments 03JAA01 & 



Table 1 ULX in the CDFN at 235.5MHz. 



Name 




s 




[mJy/beam] 


CXOHDFN J 1 2363 1 .66+620907 .3 


< 


3.63 


CXOHDFN J123632.55+621039.5 


< 


2.76 


CXOHDFN J123637. 18+621 135.0 


< 


2.88 


CXOHDFN J123641. 81+621 132.1 


< 


3.20 


CXOHDFN J123701. 47+621 845. 9 


< 


2.35 


CXOHDFN J 1 2370 1 .99+62 1122.1 


< 


3.12 


CXOHDFN J 1 23706. 1 2+62 1711.9 


< 


2.06 


CXOHDFN J123715.94+621 158.3 


< 


4.58 


CXOHDFN J123721.60+621246.8 


< 


3.47 


CXOHDFN J123723.45+621047.9 


< 


4.58 


CXOHDFN J123727.71+621034.3 


< 


4.21 


CXOHDFN J123730.60+620943.1 


< 


2.11 



ing, and tapering of the longest baselines was applied for 
N4088-X1 and N4861-X2 to improve the detection of ex- 
tended emission. 

4 Results 

4.1 ULX sources in the Chandra Deep Fields 

We show the final GMRT images obtained at 235 MHz and 
325 MHz in Fig. Q] The primary beam sizes are 1 14 ar- 
cmin at 235 MHz, 81 arcmin at 325 MHz, and 43 arcmin 
at 610 MHz. The respective rms noise in the maps at each 
frequency are 1 .4, 0.6, and 0.8 mJy/beam. No radio counter- 
parts of the ULX sources located in the CDFs are detected 
in a circle radius for each ULX position of 28 arcsec, which 
is more than 10 times the Chandra positional error circle. 
Upper limits on their flux densities at each frequency are 
given in Table 1 (for ULX located in the CDFN) and Table 
2 (ULXs in the CDFS), obtained by estimating the local rms 
at the ULX locations. These upper limits range between 2- 



1 1RNA01), observed on 2003 February 13-17 and 2007 Febru- 4 6 mJy at 235MHz, 1-2.5 mJy at 332 MHz, and 0.5-2 mJy 



ary 11-12, respectively. The observations were carried out 
using the standard GMRT phase calibration mode and spec- 
tral line mode. The data were analyzed with the NRAO As- 
tronomical Image Processing System (AIPS). In order to 
detect very faint sources, we first image the entire primary 
beam area and extract all strong point-like objects. Then the 
cleaned field of strong sources is subtracted from the visi- 
bility data, and deeper imaging is performed. Finally self- 
calibration is applied to increase the dynamic range of the 
resulting image. At 610 MHz, 3 pointings of the same field 
(CDFS) are observed. 



at 610MHz. The position of only three ULX sources fall in 
the images at 610MHz, thus only 3 upper limits are given at 
this frequency. 

4.2 Compact radio emission in ULX sources 

In Fig. |2] we show the final images of N4088-X1 (left) and 
N4861-X2 (right). The noise levels achieved are 26/iJy/beam 
for N4088-X1 and 3AtJy/beam for N4861-X2, and the restor- 
ing beams are 3 1 x 29 mas and 11x5 mas, respectively. 
For N4088-X1, we identify a compact component of 



The EVN observations (project codes EM072A & EM072BPux density 0.1 mJy at a 5er level. The component is cen- 
have been made on June 1st & 2nd 2009, in two separate tered at 



blocks of 12 hours in duration, using 9 antennas (Ef, Jb-1, 
Cm, Wb, On, Mc, Nt, Tr, Sh) at the wavelength of 18 cm. 
In order to detect weak emission from the ULX objects, we 
use the phase-referencing technique, calibrating the phases 
of the target objects with nearby strong and point-like cal- 
ibrators. The data reduction was performed in the standard 
way using AIPS. Uniform weighting was used in the imag- 



RA(J2000) = 12' l 05 m 31.7110 s ± 0.0003 s , 

DEC(J2000) = 50°32'46.729" ± 0.002". 

For this component, we estimate a brightness temperature 

ofTs>7xl0 4 K and an upper limit of 34 x 26 mas for the 

size. Adopting a distance of 13.0 Mpc yields, for N4088- 

XI, an integrated 1.6 GHz radio luminosity of 3.8 x 10 34 

erg/s. 
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GREY: HDF-N IPOL 235.500 MHZ FINAL HDFN.ICL001.1 



GREY: CDFS IPOL 332.125 MHZ FINAL CDFS.ICL001.1 
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Fig. 1 GMRT images of the HDFN at 235 MHz (left) and the CDFS at 325 MHz (right). The restoring beam sizes are 
20.1 x 16.3 arcsec at 235 MHz and 13.7 x 11.8 arcsec at 325 MHz. The best rms sensitivities achieved are 1.4, and 0.6 
mJy/beam, respectively. The ULXs positions are marked with squares. 



Table 2 

610MHz. 



ULX in the CDFS and ECDFS at 332MHz and 



Name 


S332MHZ 


SeiOMHz 




[mJy/beam] 


[mJy/beam] 


CXOECDFS J033 122.00-273620.1 


< 


1.17 




CXOECDFS J033 128.84-275904.8 


< 


1.38 


... 


CXOECDFS J033 139.05-28022 1.1 


< 


1.38 


< 1.62 


CXOECDFS J033 143.46-275527.8 


< 


1.45 




CXOECDFS J033 143.48-275 103.0 


< 


1.63 




CXOCDFS J033219. 10-274445.6 


< 


1.85 




CXOCDFS J033221.91-275427.2 


< 


2.39 




CXOCDFS J033230.01-274404.0 


< 


1.63 




CXOCDFS J033234.73-275533.8 


< 


1.28 




CXOECDFS J033249.26-273610.6 


< 


2.98 




CXOECDFS J0333 16.29-275040.7 


< 


0.99 


< 0.66 


CXOECDFS J033322.97-273430.7 


< 


1.08 


< 0.80 



The ULX N4861-X2 (Fig 12 right) has a compact com- 
ponent A centered at 

RA(J2000) = 12' l 59 m 00.3563037 ;5 ± 0.0000008 s , 
DEC(J2000) = 34°50'42.87500" ± 0.00002". 
It has a flux density of ^80/iJy (for which we derive a radio 
luminosity Li^ghz = 3.3 x 10 34 erg/s assuming a distance 
to the host galaxy of 14.80 Mpc) and a size upper limit of 
9.8 x 3.8 mas, corresponding to a brightness temperature T^ 
> 1.1 x 10 6 K. Two additional components (B and C), with 
a total flux density of ~70/iJy are detected, but cannot be 
firmly localized with the present data. If this extention were 
confirmed, the whole structure (including component A, B 



& C) would have a total flux density of 0. 18 mJy, a luminos- 
ity of Li .qghz = 7.7 x 10 erg/s and diameter D^120 mas. 
: This diameter corresponds to 8.6 pc at the distance of the 
host galaxy, and it is in agreement with the typical size of 
"HII regions found in our Galaxy, like G 18. 2-0. 3 (Fiirst et al. 
119871 ). which has a size of 200 mas, a luminosity of Li aghz 
= 1 . 1 x 10 33 erg/s and is formed by several discrete sources. 

5 Discussion 

We use the upper limits obtained from the GMRT data on 
the radio flux densities of the ULX objects to locate them 
in the fundamental plane of sub-Eddington accreting black 
holes (cf., Corbel et al. l2003l Gallo et al . 120031 Merloni et al. 
. 2003 ; Falcke et al. 2004 ) as defined by a correlation between 
radio core (£r) and X-ray (£x) luminosity and black hole 
mass, M B h, log^R = 0.6 log L x + 0.78 log M BH + 7.33. 
For our calculations, we assume a radio spectral index an ~ 
0.15 and a X-ray spectral index ax — —0.6 adopted previ- 
ously by Falcke et al. d20041 l. The resulting radio and X-ray 
luminosities of the ULX objects in our sample are compared 
in Fig.[3]to the results of Corbel et al. (2003) and Merloni et 
al. ([2003 ). The resulting high upper limits on the BH masses 
do not provide strong constraints on the nature of these ULX 
objects. 

A similar relation is shown in Fig. [4] for the most com- 
pact components in the EVN images of N4088-X1 and N4861- 
X2. Using our radio luminosity at 1.6 GHz appropriately 
scaled to 5 GHz, the X-ray luminosity scaled to the 2-10 
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CONT: N4088-X1 IPOL 1650.561 MHZ N4088-LOW.ICL001.1 
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Fig.2 1.6 GHz EVN images of the ULX sources N4088-X1 (left) and N4861-X2 (right). The restoring beam sizes are 
31 x 29 mas for N4088-X1 and 11x5 mas for N4861-X2, with the major axis of the beam oriented along a position angle 
of 21.99° for N4088-X1 and 1 1.36° for N4861-X2. The contours for the left image are (-2, 2, 3, 4, 5) x 26 ^Jy/beam, the 
rms noise off-source. For the ULX on the right, the contours are (-32, 2, 3, 4, 5, 6, 7, 10, 20, 30, 40, 50) x 3 /jjy/beam. 
These radio counterparts are offset from the X-ray peak position by 3.62" for N4088-X1, and by 1.97" for N4861-X2. 
These offsets lie within the X-ray positional error. The compact component of N4861-X2 is indicated with an A. Possible 
extended emission might be detected in regions B and C. 



keV band, and assuming a sub-Eddington accretion regime, 
we derive a black hole mass of 10 5 M and of 10 4 9 M Q 
for N4088-X1 and N4861-X2 (component A), respectively. 
These masses are in agreement with the IMBH scenario for 
both objects. 

Higher sensitivity observations are needed, and obser- 
vational time has already been guaranteed, to try to detect 
and/or confirm possible extended structure for both N4088- 
XI andN4861-X2. 

6 Summary 

Radio observations of ULX sources can help to unveil the 
nature of these objects. Analysis of archival GMRT data of 
the Chandra Deep Fields at 235, 325 & 610 MHz have not 
yielded any radio counterparts for these sources but yielded 
upper limits on their flux densities. These ULXs are too 
weak for deep field radio observations so higher sensitivity 
is needed in order to detect any faint radio emission. 

New EVN observations of three ULXs with known ra- 
dio counterparts yielded first milliarcsecond-scale images 
of all three objects. The EVN observations have confirmed 
the earlier identification of the ULX N4449-X1 with a SNR 
and obtained the most accurate estimates of its size and 
age (Mezcua & Lobanov in prep.). For the two other ULXs 
studied, N4088-X1 and N4861-X2, the EVN measurements 



have provided improved estimates of the compact radio flux, 
yielding better localizations of these objects in the Lx-L r adio 
diagram. The suggested nature of these objects can be best 
verified with more sentitive observations at 5 GHz aimed 
at both improving the brightness temperature estimates and 
obtaining spectral index information. The success of the EVN 
observations also calls for expanding this study to more ULX 
objects. 
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Fig. 3 Location of the 24 ULXs (red arrows) in the fun- 
damental plane of sub-Eddington accreting black holes. The 
parallel lines correspond to the labeled black hole mass rel- 
ative to that of the Sun. We show for comparison the Corbel 
et al. d2003l data for the X-ray binary GX 339-4 (filled cir- 
cles), and the Merloni et al. (2003) data for the Low Lumi- 
nosity AGN (LLAGN) NGC 2787, NGC 3147, NGC 3169, 
NGC 3226, and NGC 4143 (inverted triangles). 
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Fig. 4 Old (Sanchez-Sutil et al. 2006 grey squares) and 
new (this work, red squares) location of the ULX sources 
N4088-X1 and N4861-X2 in the fundamental plane of sub- 
Eddington accreting black holes. The parallel lines corre- 
spond to the labeled black hole mass relative to that of 
the Sun. We show for comparison the Corbel et al. (2003) 
data for the X-ray binary GX 339-4 (filled circles), and 
the Merloni et al. 2003 data for the Low Luminosity AGN 
(LLAGN) NGC 2787, NGC 3147, NGC 3169, NGC 3226, 
and NGC 4143 (inverted triangles). 
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